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Nepenthaceae, a monotypic family of carnivorous pitcher plants comprising Nepenthes, is widely dis¬ 
tributed in Southeast Asia. To determine the phytogeography of Nepenthes in Southeast Asia, and to 
trace the evolutionary trends of taxonomically important characteristics (i.e., peristomes) of the genus, 
we analyzed 57 internal transcribed spacer (ITS) nuclear ribosomal DNA (nrDNA) sequences of 56 spe¬ 
cies of Nepenthes and 1 ITS sequence each of Dionaea muscipula and Ancistrocladus robertsoniorum. 
To clarify the phylogenetic relationships of Nepenthes, we examined four different methods of phyloge¬ 
netic tree reconstruction. The resulting tree topologies were mostly consistent with one another except 
for the basal polytomies. Seven monophyletic subclades could be recognized. Similarities and differ¬ 
ences in terms of the positions of taxa between the present study and previous studies were observed. 
Judging from the phylogenetic trees and distribution area of each species, Borneo appears to be a second¬ 
ary center of diversification for Nepenthes and species of Nepenthes may have then radiated within the 
Sunda Shelf of Southeast Asia. The three character states of the peristomes from the upper pitchers were 
relatively well correlated with the grouping of the species of Nepenthes within seven subclades and 
showed the limitations of the Danser (1928) system for Nepenthes. 
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Nepenthaceae, a monotypic family of carniv¬ 
orous pitcher plants comprising Nepenthes, is 
widely distributed in the Asian tropics, mainly in 
Southeast Asia and the Sunda Shelf region. The 
vast majority of species grow in moist regions 
throughout the Old World tropics, as far west as 
Madagascar (TV. madagascariensis and N. maso- 
alensis ), the Seychelles (TV. pervillei), and Sri 
Lanka (TV. distillatoria ) to India (TV. khasiana ) in 
the north, Australia (TV. tenax and TV. rowanae) in 
the south; and New Caledonia (TV. vieillardii ) in 
the east (Clarke 2007, Krutzsch 1988, McPherson 
2009). 

The genus comprises 120 species in Southeast 
Asia, especially Sumatra (37 species, 29 endem¬ 
ic), Borneo (36 species, 29 endemic), and the 


Philippines (21 species, 20 endemic). Although 
most extant species of Nepenthes are distributed 
to the west of the Wallace line, some species, 
such as TV. ampullaria, TV. gracilis , TV. mirabilis, 
and TV. tentaculata, occur in both Asia and Wal- 
lacea (McPherson 2009). The distribution of the 
genus has been attributed to biogeographic fac¬ 
tors occurring both recently and in the past, in¬ 
cluding the connection of the Sunda Shelf islands 
caused by sea level drop and the isolation of the 
islands caused by increased sea levels, changes in 
global climate, and the ability of the species of 
Nepenthes to disperse and colonize new habitats 
(Clarke 2006, Danser 1928, McPherson 2009, 
Meimberg et al. 2001). 

Meimberg etal. (2001) performed a molecular 
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phylogenetic analysis of Nepenthes based on 
DNA sequences of the plastid trnK intron. The 
trnK intron is more variable than rbcL (Johnson 
& Soltis 1995, van den Berg et al. 2005) and is 
widely used for reconstruction of lower taxonom¬ 
ic levels and has therefore been the locus of choice 
in molecular phylogenetic analyses of Nepentha- 
ceae (Meimberg & Heubl 2006). The results of 
their study assumed that colonization of South¬ 
east Asia by Nepenthes was initiated from India, 
since the Indian endemic, N. khasiana, was a sis¬ 
ter taxon to all Asian taxa of Nepenthes. A com¬ 
parative analysis between the Nepenthes trnK. in¬ 
tron and its translocated copy, however, has dem¬ 
onstrated a topological incongruence. Meimberg 
et al. (2006) and Meimberg & Heubl (2006) 
pointed out the possibility of introgression or lin¬ 
eage sorting as the reason for the topological in¬ 
congruence of the phylogenetic trees. Meimberg 
& Heubl (2006) then introduced PTR1 (peptide 
transferase 1), a nuclear low copy gene as a phy¬ 
logenetic marker, but the resolution of the phylo¬ 
genetic tree was low and some taxa appeared in 
different positions from the previous findings in 
the tree topology. Meimberg et al. (2001) provid¬ 
ed two possible interpretations concerning the 
origin of Nepenthes', evolution in the northern Te- 
thys, or a Gondwanan origin at a time when the 
Indian plate was separated from Madagascar. 
From the molecular analysis based on mafK, they 
suggested that colonization of Southeast Asia 
was from an ancient Indian stock and subse¬ 
quently a new secondary center of diversity de¬ 
veloped in the Malay Archipelago. Re-evalua¬ 
tion of the origin and diversification of Nepenthes 
using another gene is needed because of topologi¬ 
cal incongruence and low resolution of the phylo¬ 
genetic trees obtained from previous studies. 

In our study, we used nucleotide sequences of 
the internal transcribed spacer (ITS) to resolve 
phylogenetic relationships within Nepenthes. 
The characteristics of the ITS region, with its 
small size, highly conserved flanking regions, 
and fast evolutionary rate, have made this nuclear 
ribosomal DNA (nrDNA) sequence a valuable 
marker for phylogenetic analysis (Baldwin et al. 
1995). In addition, ITS sequences from A. ventri- 


cosa and N. alata have been studied to evaluate 
the potential value of ITS for phylogenetic recon¬ 
struction. This study implies that the ITS regions 
of these species have many variable characteris¬ 
tics that are potentially informative for resolving 
the phylogeny of Nepenthes (Alejandro et al. 

2008) . Ribosomal DNA genes, however, are 
present in high copy numbers and may therefore 
be subjected to directional concerted evolution 
(Wendel et al. 1995a) or intergenomic introgres¬ 
sion (Wendel et al. 1995b). 

The objectives of this study were to clarify the 
phylogenetic relationships of Nepenthes based on 
ITS nucleotide sequences for 1) determining the 
phytogeography of Nepenthes in Southeast Asia 
and to re-evaluate the scenario suggested by 
Meimberg et al. (2001), and 2) trace the evolu¬ 
tionary trends of peristomes in the genus. Peri¬ 
stomes are stiff structures comprising an inward¬ 
ly curved rim surrounding the pitcher opening 
(McPherson 2009). Peristomes have been used 
in taxonomic studies (Danser 1928) to distin¬ 
guish between related species (Cheek & Jebb 
2009, Lee et al. 2006, McPherson 2009, Robinson 
et al. 2009) and to determine new species (Cheek 
& Jebb 2009, Lee at al. 2006, Robinson et al. 

2009) . 

Materials and Methods 

Plant materials 

We analyzed 57 samples from 56 species of 
Nepenthes representing all geographical areas 
(Table 1). Of the included species, seven {N. cha- 
niana, N. lingulata, N. mindanaoensis, N. naga, 
N. platychila, N. that, and N. vogelii ) were found 
within the last 10 years (Akhriadi et al. 2009, 
Cheek & Jebb 2009, Lee et al. 2006, McPherson 
2009) and have never been used in phylogenetic 
studies. Two samples from distantly distributed 
N. mirabilis were also examined. The first (N. 
mirabilisX) was from Bengkulu, Sumatra, and the 
second ( N. mirabilis2 ) was from West Kaliman¬ 
tan, Borneo. The seven species and the two dis¬ 
tantly distributed samples of N. mirabilis from 
the Malay Archipelago were expected to contrib¬ 
ute to a better understanding of the phytogeogra- 
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phy of Nepenthes in Southeast Asia. Ancistrocla- 
dus robertsoniorum (Ancistrocladaceae, Gen- 
bank: GQ443551) and Dionaea muscipula (Dro- 
seraceae, Genbank: AB675913) were used as out¬ 
groups, because these two families have been 
recognized as sister groups to Nepenthaceae on 
the basis of macromolecular characteristics using 
nuclear 18S rDNA and plastid rbcL, atpB, and 
matK sequences (Albert et al. 1992, Cuenoud et 
al. 2002, Hilu et al. 2003). 

Amplification and sequencing 

Total DNA was extracted from silica gel- 
dried leaf samples with a Qiagen DNeasy Mini 
Plant Kit (Qiagen, Valencia, CA, USA) following 
the manufacturer’s protocol. Amplification of the 
ITS region from Ancistrocladaceae was per¬ 
formed using a set of primers, AITS1 
(5’-AGAAGTCCACTGAACCTTATC-3’) and 
AITS4 (5’-CGCTTCTCCAGACTACAATTC-3’), 
which are angioseperm specific and do not co¬ 
amplify fungal DNA (Meimberg et al. 2010). The 
amplification reaction for the ITS region included 
Ex-Taq buffer and Ex-Taq DNA polymerase (Ta- 
kara Bio, Shiga, Japan). The polymerase chain 
reaction (PCR) protocol consisted of an initial 90 
s predenaturation at 96°C; 30 cycles of 20 s at 
96°C (denaturation), 40 s at 57°C (annealing), and 
40 s at 72°C (extension); and a final 7 min exten¬ 
sion at 72°C. For seven species samples (N. am- 
pullaria , N. hirsuta, N. rowanae, N. danseri, N. 
neoguineensis, N. papuana, and N. tentaculata ), 
annealing temperature was changed to 58.5°C. 

The PCR products were cleaned using the 
Wizard SV Gel and PCR Clean Up System (Pro- 
mega) and were used for autocycle sequencing re¬ 
actions following the manufacturer’s (Beckman 
Coulter) instructions. Autocycle sequencing 
products were cleaned by ethanol precipitation. 
Both forward and reverse sequences were ana¬ 
lyzed with a CEQ8000 automated sequencer 
(Beckman Coulter), using the same primers as for 
PCR. A set of internal primers, AITS2R (5’-TGC- 
GTTCAAAGACTCGATGG-3’) and AITS3F 
(5’-GAAGAACGTAGCGAAATGCG-3’), was 
designed to achieve better analysis of the ITS re¬ 
gion. The sequence of the ITS region from each 


sample was used for phylogenetic analysis; all se¬ 
quences were deposited in DDBJ/EMBL/Gen- 
bank (Table 1). 

Phylogenetic analysis 

DNA sequences obtained from the ITS region 
were aligned with ClustalX (Larkin et al. 2007). 
Phylogenetic analysis involving the maximum 
parsimony (MP) method was performed using 
the PAUP (Phylogenetic Analysis using Parsimo¬ 
ny) program, version 4.0bl0 (Swofford 2002). 
Data were analyzed by the heuristic search meth¬ 
od with the tree bisection-reconnection (TBR) 
branch swapping and MulTrees options on and 
stepwise addition with simple addition sequences 
using one reference taxon (N. that). All of the 
most parsimonious trees (MPTs) were saved. All 
characters were equally weighted and unordered 
(Fitch 1971). Gaps were treated as missing data. 
To evaluate the internal support of clades, boot¬ 
strap analysis (Felsenstein 1985) was conducted 
using 1000 replicates in a heuristic search with 
the TBR branch swapping and MulTrees options 
off. The number of steps, consistency indices, 
and retention indices (Farris 1989) were calculat¬ 
ed with one of the MPTs in each analysis using 
the TREE SCORES command in PAUP*. For 
comparison, we also performed phylogenetic 
analysis using three different methods. 

(1) Bayesian analysis with MrBayes version 
3.1.2 (Ronquist & Ehielsenbeck 2003) by using 
GTR+I+G model selected from the Modeltest 
(Posada 2008) as the optimum model for se¬ 
quence evolution based on AIC criterion. Four 
chains were run for 1,000,000 generations and 
were sampled every 100 generations to yield a 
posterior probability distribution of 10,000 trees. 
The first 2500 trees were discarded as burn in. To 
check whether the MCMC chain was long enough 
to reach convergence, the trace files resulting 
from the Bayesian analysis were opened in the 
Tracer program (vl.5) (Rambaut & Drummond 
2007), and their effective sample size (ESS) sta¬ 
tistics were calculated for values higher than 100. 

(2) The neighbor joining (NJ) method (Saitou 
& Nei 1987) with MEGA version 5.05 (Tamura et 
al. 2011) by using the Tamura 3-parameter model 
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Table 1. Plant materials examined in this study. The materials were collected from the Indonesian Carnivorous Plant Society 
(Komunitas Tanaman Karnivora Indonesia/KTKI), Bogor Botanical Garden (BBG, Indonesia), and Kyoto Botanical Gar¬ 
dens (KBG, Japan). All specimens were deposited in TI (*FA: Firman Alamsyah). 

Species Source Voucher Voucher No. of Genbank 

_ Specimens* living stocks accession no. of ITS 

Nepenthes adnata Tamin & M. Hotta ex Schlauer 


Nepenthes alata Blanco 
Nepenthes alba Ridl. 

Nepenthes ampullaria Jack 
Nepenthes bellii K. Kondo 
Nepenthes burbidgeae Hook. f. ex Burb 
Nepenthes burkei Mast. 

Nepenthes campanulata Sh. Kurata 
Nepenthes chaniana C. Clarke et al. 

Nepenthes clipeata Danser 
Nepenthes copelandii Merr. ex Macfarlane 
Nepenthes danseri Jebb & Cheek 
Nepenthes densiflora Danser 
Nepenthes diatas Jebb & Cheek 
Nepenthes distillatoria L. 

Nepenthes ephippiata Danser 
Nepenthes faizaliana J.H. Adam & Wilcock 
Nepenthes fusca Danser 

Nepenthes glabrata J.R. Turnbull & A.T. Middleton 
Nepenthes gracilis Korth. 

Nepenthes hirsuta Hook. f. 

Nepenthes khasiana Hook. f. 

Nepenthes lingulata Chi. C. Lee, Hernawati & Akhriadi 
Nepenthes longifolia Nerz & Wistuba 
Nepenthes macrovulgaris J.R. Turnbull & A.T. Middleton 
Nepenthes madagascariensis Poir. 

Nepenthes merrilliana Macfarlane 
Nepenthes mindanaoensis Sh. Kurata 
Nepenthes mirabilis (Lour.) Druce (Bengkulu, Sumatra) 
Nepenthes mirabilis (Lour.) Druce (West Kalimantan) 
Nepenthes naga Akhriadi et al. 

Nepenthes neoguineensis Macfarlane 
Nepenthes ovata Nerz Sc Wistuba 
Nepenthes papuana Danser 
Nepenthes pervillei Blume 
Nepenthes platychila Chi. C. Lee 
Nepenthes rajah Hook. f. 

Nepenthes reinwardtiana Miq. 

Nepenthes rhombicaulis Sh. Kurata 
Nepenthes rowanae F.M. Bailey 
Nepenthes sanguinea Lindl. 

Nepenthes smilesii Hemsl. 

Nepenthes spathulata Danser 
Nepenthes spectabilis Danser 
Nepenthes sp. Misool 
Nepenthes stenophylla Mast. 

Nepenthes sumatrana (Miq.) Beck 
Nepenthes talangensis Nerz & Wistuba 
Nepenthes tentaculata Hook. f. 

Nepenthes thai Cheek 
Nepenthes tobaica Danser 
Nepenthes truncata Macfarlane 
Nepenthes veitchii Hook. f. 

Nepenthes ventricosa Blanco 
Nepenthes vieillardii Hook. f. 

Nepenthes villosa Hook. f. 

Nepenthes vogelii Schuit. & de Vogel 
Dionaea muscipula J. Ellis ex L. 


KTKI 

FA-110001 

B111109 01 

AB675864 

KTKI 

FA-110002 

C14110901 

AB675865 

KTKI 

FA-110003 

E05021001 

AB675866 

BBG 

FA-110049 

KRB27100901 

AB675914 

KTKI 

FA-110004 

Al 11109 01 

AB675868 

KTKI 

FA-110005 

B14110901 

AB675869 

KTKI 

FA-110006 

E05021002 

AB675870 

KTKI 

FA-110007 

B14110902 

AB675871 

KTKI 

FA-110008 

Bill 10902 

AB675872 

KTKI 

FA-110009 

D07110901 

AB675873 

KTKI 

FA-110010 

E05021003 

AB675874 

KTKI 

FA-110050 

A14110901 

AB675915 

KTKI 

FA-110057 

C14110902 

AB675875 

KTKI 

FA-110011 

D071109 02 

AB675876 

KTKI 

FA-110012 

E05021004 

AB675877 

KTKI 

FA-110013 

A14110902 

AB675878 

KTKI 

FA-110014 

E05021005 

AB675879 

KTKI 

FA-110015 

A141109 03 

AB675880 

KTKI 

FA-110016 

All 110902 

AB675881 

BBG 

FA-110017 

KRB27100902 

AB675882 

KTKI 

FA-110051 

D07110903 

AB675916 

KTKI 

FA-110018 

E05021006 

AB675883 

KTKI 

FA-110019 

Al 1110903 

AB675884 

KTKI 

FA-110020 

E05021007 

AB675885 

KTKI 

FA-110021 

B14110903 

AB675886 

KBG 

FA-110058 

KBG02-0521 

AB769064 

KTKI 

FA-110022 

A11110904 

AB675887 

KTKI 

FA-110023 

A141109 04 

AB675888 

KTKI 

FA-110024 

E05021008 

AB675889 

KTKI 

FA-110025 

E05021009 

AB675890 

KTKI 

FA-110026 

D07110904 

AB675891 

KTKI 

FA-110052 

A14110905 

AB675917 

KTKI 

FA-110027 

Bill10903 

AB675892 

KTKI 

FA-110053 

B141109 04 

AB675918 

KTKI 

FA-110028 

E05021010 

AB675893 

KTKI 

FA-110029 

E05021011 

AB675894 

KTKI 

FA-110030 

D07110905 

AB675895 

BBG 

FA-110031 

KRB27100903 

AB675896 

KTKI 

FA-110032 

B141109 05 

AB675897 

KTKI 

FA-110054 

A14110906 

AB675919 

KTKI 

FA-110033 

B14110906 

AB675898 

KTKI 

FA-110034 

C14110903 

AB675899 

KTKI 

FA-110035 

D07110906 

AB675900 

KTKI 

FA-110036 

A14110907 

AB675901 

KTKI 

FA-110037 

E05021012 

AB675902 

KTKI 

FA-110038 

A141109 08 

AB675903 

KTKI 

FA-110039 

A14110909 

AB675904 

KTKI 

FA-110040 

D07110907 

AB675905 

KTKI 

FA-110055 

C141109 04 

AB675920 

KTKI 

FA-110041 

E05021013 

AB675906 

KTKI 

FA-110042 

Bill10904 

AB675907 

KTKI 

FA-110043 

E05021014 

AB675908 

KTKI 

FA-110044 

A14110910 

AB675909 

KTKI 

FA-110045 

Al 1110905 

AB675910 

KBG 

FA-110059 

KBG02-0522 

AB769065 

KTKI 

FA-110046 

B14110907 

AB675911 

KTKI 

FA-110047 

B11110905 

AB675912 

KTKI 

FA-110048 

E05021015 

AB675913 
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(Tamura 1992). Bootstrap values were calculated 
with 1000 replicates. 

(3) The maximum likelihood (ML) method 
with the Treefinder program (Jobb et al. 2004) by 
using GTR+G model, which was suggested based 
on AIC criterion by the “propose model” analysis 
implemented in the Treefinder program. Boot¬ 
strap values were calculated with 1000 replicates. 

Reconstruction of character states 

To study the phytogeography of Nepenthes in 
Southeast Asia, we mapped the character states 
for distribution areas onto one of the MPTs by us¬ 
ing the MacClade program (v4.06), with acceler¬ 
ated character transformation (ACCTRAN) opti¬ 
mization (Maddison & Maddison 2003). In addi¬ 
tion, we also mapped the character states for the 
peristome of upper pitchers to study the evolu¬ 
tionary trends of morphological characteristics of 
Nepenthes. The morphometric data of the peri¬ 
stomes of the upper pitchers were taken from 
McPherson (2009). 

Upper pitchers are distinguishable from lower 
pitchers. They are produced from leaves along a 
climbing stem of older plants of Nepenthes, in 
which the tendril attaches to the back of the pitch¬ 
er on the side where the lid meets the rear of the 
pitcher orifice. In contrast, the lower pitchers, 
which are developed from a tendril that attaches 
to the front of the pitcher on the same side as the 
orifice, are produced by young plants of Nepen¬ 
thes (Di Giusto et al. 2008, Gaume & Di Giusto 
2009, McPherson 2009, Moran 1996). All but 
one of the species of Nepenthes examined in this 
study showed this pitcher dimorphism. Nepen¬ 
thes campanulata produces only one type of 
pitcher. We included data for the peristomes of 
these pitchers, since they develop from leaves 
farther up the stem of mature plants and show 
tendril attachment as in the typical upper pitch¬ 
ers of Nepenthes. We classified peristomes as 
having three character states: narrow (<10 mm), 
intermediate (10-20 mm), and broad (>20 mm). 

Results 

Molecular phylogeny o/Nepenthes 


We obtained 57 new ITS nrDNA sequences of 
56 species of Nepenthes and 1 ITS sequence of 
Dionaea muscipula. The alignment of 59 entire 
ITS sequences (57 ingroup sequences plus 2 out¬ 
groups) provided an 856-bp-long matrix. Se¬ 
quence length variations resulting from inser¬ 
tions and deletions were found among the species 
of Nepenthes. The aligned ITSs contained 341 
(39%) constant characters, 239 (29%) parsimony- 
uninformative variable characters, and 276 (32%) 
parsimony-informative characters. The analysis 
resulted in 324 MPTs with a length of 1130 steps 
and had consistency (Cl) and retention (RI) indi¬ 
ces of 0.665 and 0.687, respectively. The strict 
consensus tree reconstructed by the parsimony 
method is shown in Fig. 1. The trees obtained 
from the NJ, ML, and Bayesian analysis methods 
were mostly consistent with the tree obtained 
from the parsimony method, except for basal spe¬ 
cies, which appeared in different positions (Fig. 1 
& Fig. 2). We compiled the bootstrap values of 
MP, NJ, ML, and posterior probabilities of Bayes¬ 
ian analysis on the MP tree (Fig. 1). 

Using Ancistrocladus robertsoniorum and 
Dionaea muscipula as outgroups, we recognized 
two basal branches and seven subclades (I-VII) 
in the phylogenetic tree (Fig. 1). In addition, we 
recognized a basal group consisting of seven spe¬ 
cies whose branching pattern could not be re¬ 
solved well. The seven subclades (I-VII) were 
supported by all four phylogeny reconstruction 
methods, except for subclade I. The members of 
subclade I were grouped with some of the basal 
unresolved species in the Bayesian tree (Fig. 2). 
Corresponding to the tree topology, the two basal 
taxa were N. pervillei from Seychelles and N. 
madagascariensis from Madagascar. Subclade I 
comprised two species from outlying areas (India 
(N. khasiana ) and Sri Lanka (N. distillatoria)); a 
species from New Guinea (N. papuana ), and a 
species (N. ampullaria ) distributed across Penin¬ 
sular Malaysia, Sumatra, Borneo, and New Guin¬ 
ea. Subclade II consisted of a species endemic to 
Sulawesi (N. glabrata), a species from Sulawesi 
and Borneo (N. tentaculata ), and a species en¬ 
demic to Borneo (N. hirsuta). Subclade III con¬ 
tained species restricted to the Philippines and a 
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species endemic to Borneo (N. campanulata ) in 
the basal position. Subclade IV also contained 
species restricted to the Philippines. Subclade V 
contained only species from Borneo. Subclade 
VI comprised exclusively four species from Pen¬ 
insular Malaysia and Indochina (N. smilesii, N. 


sanguinea, N. alba, and N. thai ), a species en¬ 
demic to Borneo (N. macrovulgaris ), and a spe¬ 
cies occuring in Borneo and Sumatra ( N . rein- 
wardtiana) in the basal position. Subclade VII 
comprised 13 species of Nepenthes from Sumatra 
and one species ( N . spathulata ) also distributed in 



Sey: Seychelles 

Mad: Madagascar 

Ind: India 

Sri: Sri Lanka 

Aus: Australia 

NC: New Caledonia 

NG: New Guinea 

Sul: Sulawesi 

Bor: Borneo 

Phi: Philippines 

MI: Malaysia and Indochina 

Sum: Sumatera 


71/61/1 
67/0.98 87/57*' 
/ 86 / 1.00 


muscipula 
robertsoniorum 
pervillei (Sey) 
madagascamnsis (Mad) 
neoguineensis (NG) 
danseri (NG) 
rowanae (Aus) 
vieillardii (NC) 
gracilis (Sul,Bor,Ml,Sum) 
sp. Misool (NG) 

mirabilisl (Aus,NG,Sul,Bor,Phi.MI,Sum) 

mirabilis2 (Aus.NG.SufBor.Phi.MfSum) 

distillatoria (Sri) 

khasiana (ind) 

papuana (NG) 

ampuilaria (NG,Bor,MI,Sum) 

glabrata (Sul) 

tentaculata (Sul,Bor) 

hirsuta (Bor) 

campanulata (Bor) 

mindanaoensis (Phi) 

copelandii (Phi) 

truncata (Phi) 

alata (Phi) 

beilii (Phi) 

merrilliana (Phi) 

ventricosa (Phi) 

burkei (Phi) 

villosa (Bor) 

rajah (Bor) 

veitchii (Bor) 

dipeata (Bor) 

stenophytla (Bor) 

chaniana (Bof) 

ephippiata (Bor) 

faizaliana (Bor) 

fusca (Bor) 

platychila (Bor) 

burbidgeae (Bor) 

vogelii (Bor) 

reinwardtiana (Bor,Sum) 
macrovulgaris (Bor) 
thai (Ml) 
sanguine a (Mi) 
alba (Ml) 
smilesii (Ml) 
longifolia (Sum) 
naga (Sum) 
adnata (Sum) 
rhombicauiis (Sum) 
lingulata (Sum) 
ova la (Sum) 
talangensis (Sum) 
spathulata (Sum) 
sumatrana (Sum) 
tobaica (Sum) 
densitlora (Sum) 
diatas (Sum) 
spectabilis (Sum) 


Basal 

species 


III 


IV 


V 


Vi 


Vli 


Fig. 1. Strict consensus tree derived from maximum parsimony analysis of ITS sequences of Nepenthes and outgroup taxa. 
Statistical support for each branch is shown below each branch with successive values of MP/NJ/ML/Bayesian. The mi¬ 
nus signs (-) on the successive statistical values indicate different topologies or polytomies. Additional information on the 
distribution area of a the species is given as an abbreviation next to the taxon name. 
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N. adnata 
N. rhombicaulis 
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N. ovata 
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N. diatas 
N. spectabilis 
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Fig. 2. The 50% majority-rule consensus tree derived from the Bayesian analysis of ITS sequences for Nepenthes and outgroup 
taxa. Numbers below branches are Bayesian posterior probabilities. 
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Fig. 3. Character state reconstruction of Nepenthes for distribution areas based on 1 of 324 MPTs using MacClade ver. 4.06 
with ACCTRAN optimization. The squares next to the taxon name indicated subclades I—VII. 
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Fig. 4. Character state reconstruction of Nepenthes for peristomes based on 1 of 324 MPTs using MacClade ver. 4.06 with 
ACCTRAN optimization with comparison to the classification system of Nepenthes by Danser (1928), which is indicated 
by the squares containing group names next to the taxon names. The arrows indicate repeated evolution of narrow, inter¬ 
mediate, and broad peristomes. We recognized three character states of the peristomes: narrow (<10 mm), intermediate 
(10-20 mm), and broad (>20 mm). The character state for each species was obtained from the descriptions by McPherson 
(2009). 
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Java. 

Seven species were indicated as basal polyto- 
mies on the MP phylogenetic tree, including two 
nonendemic species (N. mirabilis and N. gracilis) 
distributed across Borneo, Sumatra, Sulawesi, 
Peninsular Malaysia and Indochina and the Phil¬ 
ippines, New Guinea, and northern Australia for 
N. mirabilis', a species from Misool Island located 
near the west coast of New Guinea (N. sp. Misool), 
two species from New Caledonia (N. vieillardii ) 
and Australia (N. rowanae), and two species from 
New Guinea ( N. neoguineensis and N. danseri), 
and Waigeo Island located near the northwest 
coast of New Guinea (N. danseri). 

Character state reconstruction 

The character state reconstructions of Nepen¬ 
thes for distribution areas and peristomes, based 
on one of the 324 MPTs are presented in Fig. 3 
and Fig. 4, respectively. Figure 3 shows the evo¬ 
lutionary trends of distribution areas of Nepen¬ 
thes, which exhibits the radiation of Nepenthes on 
some islands and the migration of Nepenthes to 
adjacent islands or to mainland of Southeast Asia. 
Figure 4 shows the evolutionary trends of the 
peristomes of the upper pitchers and a compari¬ 
son of the characteristics of the peristomes with 
the classification system of the genus by Danser 
(1928). 

Discussions 

Phylogenetic relationships in the genus Nepen¬ 
thes 

In our study of ITS DNA analysis (Fig. 1), we 
obtained phylogenetic trees with relatively higher 
bootstrap supports and Bayesian posterior prob¬ 
abilities than in previous studies employing the 
DNA region in the chloroplast genome (Meim- 
berg et al. 2001) or the coding region of the nucle¬ 
ar genome (Meimberg & Heubl 2006), except for 
the basal positions of the trees. It is advantageous 
to use nuclear ITS regions for phylogenetic analy¬ 
sis in angiosperms. 

Nepenthes pervillei was found to be the most 
basal taxon within the genus. It can be distin¬ 
guished from all other species by its seeds, which 


lack the appendages typical of most species of 
Nepenthes. In addition, N. pervillei also has 
some unusual characteristics not found in other 
species of Nepenthes, such as nontwining tendrils 
that only slightly emerge from the upper pitcher; 
similar upper and lower pitchers; black, short, 
ovoid, truncate seeds and obconic fruit (Meim¬ 
berg et al. 2001, McPherson 2009). The basal po¬ 
sition of N. pervillei was consistent with reports 
from previous studies of the phylogeny of Nepen¬ 
thes using the trnK intron (Meimberg et al. 2001, 
Meimberg & Heubl 2006, Meimberg et al. 2006). 
The second basal taxon in our study, N. madagas- 
cariensis, was also consistent with the trnK. in¬ 
tron phylogeny of Nepenthes (Meimberg et al. 
2001, Meimberg et al. 2006), but in a subsequent 
study of the phylogeny of Nepenthes using PTR1, 
N. madagascariensis formed a clade together 
with an Indian species, N. khasiana, and a Su- 
lawesian species, N. tomoriana (Meimberg & 
Heubl 2006). 

In our study, all the New Guinean species 
were included with the basal species and in sub- 
clade I with species from outlying areas, includ¬ 
ing Australia (N. rowanae), India (N. khasiana), 
and Sri Lanka (TV. distillatoria), and the widely 
distributed N. ampullaria. The positions of the 
New Guinean species N. papuana in the same 
subclade as species from India (TV. khasiana ) and 
Sri Lanka (TV. distillatoria) differ from previous 
studies of the trnK intron phylogeny of Nepen¬ 
thes, where the Indian species ( TV . khasiana) was 
recognized as the third basal branch. The posi¬ 
tions of the widely distributed N. gracilis and N. 
mirabilis as basal species, and N. ampullaria in 
subclade I of the tree topology as determined in 
our study were also different from reports in pre¬ 
vious studies of the trnK intron phylogeny of Ne¬ 
penthes, where these three widely distributed 
species were included in the same subclade as the 
Bornean species (Meimberg et al. 2001, Meim¬ 
berg & Heubl 2006, Meimberg et al. 2006). 

Despite the different positions of some spe¬ 
cies in subclade I and the basal species with pre¬ 
vious studies of Nepenthes phylogeny, these spe¬ 
cies share similar morphological characteristics 
including acute leaf apex, orbicular or partly or- 
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bicular pitcher lids, and partly cylindrical lower 
pitcher form, except for N. ampullaria, which has 
a distinctive urceolate shape (McPherson 2009). 
Moreover, the pitcher morphology of N. mirabilis 
was shared by three species from outlying areas 
(N. khasiana , N. distillatoria, and N. vieillardii), 
which similarly lack putative apomorphic charac¬ 
teristics (McPherson 2009). 

The positional differences of species between 
the present ITS analysis and previous studies us¬ 
ing P-nK intron and PTR1 (Meimberg et al. 2001, 
Meimberg & Heubl 2006, Meimberg et al. 2006) 
may be due to differences in the markers used. 
ITS is located in the nuclear genome (Baldwin et 
al. 1995), whereas the trnK intron is located in 
the chloroplasts, which are maternally inherited 
and may lead to the chloroplast capture phenom¬ 
enon (Soltis & Kuzoff 1995). This phenomenon 
is considered the main reason for the incongru¬ 
ence of nuclear markers and plastid phylogenies 
(Meimberg & Heubl 2006, Tsitrone et al. 2003). 
In contrast, ITS has another weaknesses for esti¬ 
mating phylogeny. ITS is included in the rRNA 
precursor transcript. Genes encoding rRNA and 
spacers occur in tandem repeats that are thou¬ 
sands of copies long and polymorphisms in their 
sequences are often observed (Wendel et al. 
1995a). Indeed, in the species of Nepenthes ex¬ 
amined, we detected polymorphisms in seven 
species when using the PCR protocol with an¬ 
nealing temparature at 57°C. It is critical to de¬ 
cide which is an orthologous sequence for recon¬ 
structing the phylogenetic trees. It is possible 
that the incongruence of phylogenetic tree topol¬ 
ogies is due to using paralogous ones. The con¬ 
certed evolution (Wendel et al. 1995a) and in- 
tergenomic introgression (Wendel etal. 1995b) of 
the nrDNA genes may also have affected this in¬ 
congruence. In the case of PTR1 (Meimberg & 
Heubl 2006), the gene is located in the nuclear ge¬ 
nome, the same as the ITS. The incongruence 
between ITS and PTR1 trees may be due to am¬ 
plification of paralogous sequences. As men¬ 
tioned above, ITS sometimes shows polymor¬ 
phisms, and PTR1 also exhibits the possibility of 
paralogous genes, since Nepenthes is suspected 
to be of polyploid origin (Meimberg & Heubl 
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2006). Lineage sorting is also considered to be a 
candidate for this incongruence. In addition to 
the incongruence discussed above, the basal re¬ 
gion of our trees was not resolved well, and it is 
highly desirable to examine the phylogeny of Ne¬ 
penthes using other DNA regions. 

Phytogeography o/Nepenthes in Southeast Asia 

Meimberg et al. (2001) suggested the scenario 
that Nepenthes colonized Southeast Asia from an 
ancient Indian stock, first to the Malay Peninsula, 
and subsequently to Indochina and the Malay Ar¬ 
chipelago. Their results showed the Indian spe¬ 
cies to be sister to the Southeast Asian species, 
although they did not reach a final conclusion. 
Our results show the Indian and Sri Lanka spe¬ 
cies forming a clade with the New Guinean N. 
papuana and with the widely distributed N. am¬ 
pullaria and does not support the suggested sce¬ 
nario of Meimberg et al. (2001). In addition, three 
other New Guinean species, N. sp. Misool, N. 
neoguineensis, and N. danseri, were included in 
the basal polytomies with other widely distribut¬ 
ed species, i.e., N. gracilis and N. mirabilis. The 
distribution area reconstruction on one of the 324 
MPTs (Fig. 3) shows that these three widely dis¬ 
tributed species may have expanded from east to 
west in Southeast Asia (east to west hypothesis). 
Of course, due to the polytomy of the basal spe¬ 
cies, the scenario of Meimberg et al. (2001) is 
also possible based on other MPTs. They consid¬ 
ered that expansion of Nepenthes occurred from 
west to east in Southeast Asia, based on their tree 
topology, in which Indian species N. khasiana 
was sister to all of the Southeast Asian species. 
That topology was not supported by the PTR1 
tree, where N. khasiana formed a clade with the 
species of Sulawesi and Madagascar (Meimberg 
& Heubl 2006). The east to west hypothesis was 
corroborated by similarities in morphology, alti¬ 
tudinal distribution and habitat types of the three 
widely distributed species, which were shared by 
species from New Guinea, Australia and New 
Caledonia (McPherson 2009). 

The two species from Sulawesi, N. glabrata 
and N. tentaculata were included in subclade II 
together with N. hirsuta, a species endemic to 
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Borneo. This subclade is sister to subclade III" 
VII consisting of most species occurring in Bor¬ 
neo, Sumatra, the Philippines, Malaysia and In¬ 
dochina. This tree topology and presences of N. 
tentaculata in Sulawesi and the Bornean Islands 
also suggests an important linkage for species of 
Nepenthes distributed across Southeast Asia 
from east to west (Fig. 3). The flora of Sulawesi is 
closely linked with the Australian-New Guinean 
phytogeographic region, which is indicated by 
the floral exchange between these three areas 
(van Welzen etal. 2011). Therefore, Sulawesi may 
have been a bridge for the migration of Nepenthes 
from east to west. 

The distribution area reconstruction of Ne¬ 
penthes (Fig. 3) suggests that Nepenthes diverged 
in Borneo, resulting many local species and could 
have also migrated to adjacent areas. Nepenthes 
probably migrated from Borneo to the Philip¬ 
pines. This scenario is supported by the presence 
of N. campanulata in the basal position of sub¬ 
clade III, which includes Philippine species. In 
addition, subclade IV consisting of Philippines 
species was included in clades (subclade III—VII), 
which is sister to subclade II consisting of species 
from Borneo and Sulawesi. This scenario also 
differs from that of Meimberg et at. (2001), which 
suggested migration from the Philippines to Bor¬ 
neo. In subclade VI, the two Bornean species 
were located in basal positions as sisters to the 
Indochina species and to all species of subclade 
VII that occur on Sumatra. Our results therefore 
suggest an ancestral species migrated from Bor¬ 
neo, diverged in Sumatra, then subsequently mi¬ 
grated to Java (N. spathulata is also known to be 
distributed in Java). It can be assumed that Bor¬ 
neo was a secondary center of diversification for 
Nepenthes , allowing species of Nepenthes to ra¬ 
diate within the Sunda Shelf region. 

The scenario for Nepenthes expansion esti¬ 
mated from our molecular analysis based on ITS 
differs from that of Meimberg et al. (2001). 
These two scenarios were based on results from 
different DNA regions located in different ge¬ 
nomes and may be due to the characteristics of 
the different DNA regions used in the phyloge¬ 
netic analyses. The species examined in both 
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analyses also differed, possibly affecting the tree 
topologies. Because the phylogenetic analysis 
based on ITS has weaknesses, especially in dis¬ 
criminating between orthologs and paralogs, it is 
highly desirable to base phylogenetic reconstruc¬ 
tions on other DNA regions and by using addi¬ 
tional species of Nepenthes. 

Morphology of the peristome 

Some morphological characteristics in the 
species of Nepenthes are polymorphic, and each 
of them has more than one form or shape within a 
species, making it difficult to evaluate the evolu¬ 
tionary trends in morphological characteristics. 
In this study, we tried to analyze some morpho¬ 
logical characteristics, such as indumentum, in¬ 
florescence, lamina, leaf apex and base, lower 
and upper pitcher forms and the lid and peristome 
of the pitcher. We found the peristome of the up¬ 
per pitchers to be consistent, except in N. mindan- 
aoensis, N. alata, N. ovata, and N. densiflora. 
The peristome of the upper pitchers were relative¬ 
ly well correlated with the grouping of Nepenthes 
into seven subclades. 

Figure 4 shows the evolutionary trends in the 
peristomes. The peristomes of the two basal spe¬ 
cies, as well as in most species in the basal poly- 
tomies and subclades I, II, and III, were narrow, 
suggesting that the narrow peristome is the ple- 
siomorphic state in Nepenthes. In addition, all 
species in subclade VI produced narrow peri¬ 
stomes. Whereas the Bornean species in sub¬ 
clade V and most species in subclade IV and VII, 
as well as N. ampullaria in subclade I and two 
basal species (N. neoguineensis and N. rowanae ), 
had intermediate and broad peristomes. The in¬ 
termediate and broad peristomes likely evolved at 
least seven or eight times, respectively, whereas 
narrow peristomes evolved at least five times af¬ 
ter the widening of the peristome. The findings 
suggest repeated evolution of narrow, intermedi¬ 
ate, and broad peristomes. 

In the taxonomic study by Danser (1928), nar¬ 
row peristomes were considered to be ancestral 
in the Vulgatae group, which is concordant with 
our results, suggesting that N. pervillei and N. 
madagascariensis are basal species and that the 
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tree topology of subclade I and the basal polyto- 
mies consist of three species from the outlying 
areas (N. distillatoria, N. khasiana, and N. vieil- 
lardii ), two New Guinean species (N. danseri and 
N. papuana ), and two widely distributed species 
(N. gracilis and N. mirabilis). All of these species 
were assigned by Danser (1928) to the Vulgatae 
group (Fig. 4). Although N. neoguineensis in the 
basal polytomies produces intermediate peri¬ 
stomes, Danser (1928) assigned it to the Vulgatae 
group. Other species with narrow peristomes 
from three different subclades were also assigned 
by Danser (1928) to the Vulgatae group. Since 
representatives of the Vulgatae group are scat¬ 
tered in different subclades (Fig. 4), the group 
with narrow peristomes is obviously polyphylet- 
ic. 

Intermediate peristomes can be defined as 
one of the key characteristica of the Nobiles 
group, while broad peristomes can be defined as 
a characteristic of the Regiae, Insignes and Ur- 
ceolatae groups (Danser 1928). Since represen¬ 
tatives of the Insignes, Regiae, and Nobiles 
groups occur in different subclades (Fig. 4), the 
three groups that share intermediate and broad 
peristomes are obviously polyphyletic. Judging 
from the tree topology (Fig. 4), the evolutionary 
trends in the character states of the peristomes re¬ 
veal the limitations of Danser’s classification of 
Nepenthes. 
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